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Abstract 

Flight measurements of pitch, yaw, and roll rates and the resultant 
rotational ly Induced'nnear accelerations during three orbital maneuvers on 
Shuttle mission STS 61-C have been used to calculate the actual orbiter 
center-of-gravity location. The calculation technique reduces error 
due to lack of absolute calibration of the accelerometer measurements and 
compensates for accelerometer temperature bias and for the effects of 
gravity gradient. Accuracy of the technique was found to be limited by the 
nonrandom and asymmetrical distribution of orbiter structural vibration at 
the accelerometer mounting location. Fourier analysis of the vibration was 
performed to obtain the power spectral density profiles which show 
magnitudes in excess of 10 4 ug 2 /Hz for the actual vibration and over 
500 ug 2 /Hz for the filtered accelerometer measurements. The data from this 
analysis provide a characterization of the Shuttle acceleration environment 
which may be useful in future studies related to accelerometer system 
application and "zero-g" investigations or processes. 
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g 
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Mg 

PSD 


P, Q, r 


P, Q, r 


coefficients for HiRAP temperature bias correction 
(i = x,z) 

acceleration induced by gravity gradient (i = x,z) 
total Induced acceleration (i = x,z) 

measured acceleration (i = x,z) 

acceleration induced by rotation (i = x,y,z) 
correction for zero offset and temperature bias (i = x, 

center of gravity 
Greenwich Mean Time 

acceleration of gravity at sea level (9.806 m/sec^) 
orbiter moments of inertia about X^, Y 5 , Zb 
orbiter pitching moment induced by gravity gradient 
power spectral density 

angular velocities about orbiter Xb, Yb» Zb axes 
angular acceleration about orbiter Xb, Yb, Zb axes 


q 9 

R 

t 

tref 

tl.t 2 

t 3 ,t 4 
X, Y,Z 
X b ,Y b ,Z b 


pitch angular acceleration induced by gravity gradient 

distance from Earth to orbiter 

time 

reference time for temperature bias correction 
start and end times for pre-manuever data segment 
start and end times for manuever data segment 
HiRAP axes (parallel to orbiter body axes) 
orbiter body axes 
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Xi ,Yi ,Zi 

e 

o 

w 


offsets of the HiRAP X,Y,Z accelerometer from the 
orbiter center of gravity along the ith body axis 
(i = x,y,z) 

preflight values of X-j, Yj, and Zi 

angle between X^ and local vertical 

Earth gravitational constant (3.989xl0 14 m 3 /sec 2 ) 

standard deviation 

orbital angular velocity 


Introduction 

Measurement of very low level linear acceleration of the Shuttle Orbiter 
is one of the principal experiments being conducted as part of the Orbiter 
Experiments (OEX) Project. The High Resolution Accelerometer Package (HiRAP) 1 
has been providing flight measurements since the sixth Shuttle flight on all of 
the 0V-102 Columbia and OV-99 Challenger flights. HiRAP is capable of 
measurement in the low 10 - ® g range, which has facilitated the determination of 
orbiter aerodynamic coefficients and atmospheric density in the hypersonic, 
rarefied-flow transition region. An even more sensitive accelerometer system, 
the Orbital Acceleration Research Experiment (OARE) 2 , is ready for Installation 
and will return acceleration measurements in the low lO - ^ g range. Accurate 
interpretation of the OARE flight data will require very accurate knowledge of 
the actual orbiter center-of-gravi ty location at the time of measurement. This 
requirement is the stimulus for this analysis of the possibility of orbiter 
center-of-gravity determination from flight data. Additionally, the technique 
herein or the results from this study may be of interest relative to attitude 
guidance and control of large-scale space structures or to experiments 
sensitive to the orbiter "gravity" environment. 
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This analysis is formulated on the concept of induced acceleration sensed 
by a linear accelerometer when the body to which the accelerometer Is mounted 
rotates about its center of gravity. The induced acceleration is proportional 
to the angular velocity of the body and the offset of the accelerometer axis 
from the center of gravity. The flight measurements of acceleration for the 
analysis were obtained from the HiRAP, and the angular velocity data were 
obtained from the rate gyros on the Aerodynamic Coefficient Identification 
Package (ACIP), a companion OEX experiment. This report discusses the 
dedicated orbiter attitude maneuvers for the experiment, the resultant flight 
data, development of the analysis technique, and results of the analysis. 

STS 61-C OEX Attitude Maneuvers 

Three attitude maneuvers were performed during the orbital phase of the 
STS 61-C mission flown by the 0V-102 Orbiter (Columbia) in January 1986. These 
maneuvers were designed to support several OEX experiments. Including this 
analysis of in-flight center-of-gravity determination. The maneuver 
specifications 3 were intended to produce a simple pitch rotation with minimum 
forces and moments during execution of the maneuver. Figure 1 shows the 
maneuver geometry with the orbiter initially belly forward and nose down at an 
angle of 20° with the local vertical. This pre-maneuver attitude was 
maintained for several minutes by holding a positive pitch rate equal to the 
orbital angular velocity of 0.0658 deg/sec. The maneuver was initiated with 
the firing of coarse attitude thrusters to achieve a negative pitch rate of 
about 0.5 deg/sec, after which the orbiter was allowed to coast in the drift 
mode (no thruster firings) to the final position with belly aft and nose up at 
180° from local vertical. Yaw and roll rates were held to very small values 
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with ±5° tolerances on the target values of 0° for yaw and roll angles. Crew 
motion, dumps, and vents were inhibited during the maneuvers. ACIP and HiRAP 
data were recorded for several minutes during the pre-maneuver attitude hold, 
during the pitch rotation, and for several minutes during the post-maneuver 
attitude hold. 


ACIP Rate Gyro Flight Data 

The flight data from the ACIP rate gyros were recorded during the attitude 
maneuvers and processed postflight at the OEX Data Lab (ODL) at NASA Johnson 
Space Center (JSC). The processing included calibration corrections for zero 
bias, temperature, cross-coupling, and other error sources. ^ The final 
calibrated data were provided to NASA Langley Research Center (LaRC) in the 
standard OEX computer-compatible tape (CCT) format at a data sample rate of 
112.7 per second. The pitch rate measured by the q-channel rate gyro is shown 
on figures 2-4 for the three maneuvers and generally conforms to the target 
specifications with only three noticeable exceptions. The initial overshoot on 
the pitch rate at the start of maneuver 1 was followed by a very slow 
adjustment to acquire the final -0.67 deg/sec rate. This adjustment appears to 
have been obtained by pulsing the vernier thrusters over a 30-second interval. 

A second deviation occurred during maneuver 3 which shows an initial rate 
during the pre-maneuver attitude hold which was below the orbital rate, 
resulting in a slow drift toward local vertical. A correction was made with 
the vernier thrusters at about t = 55640 seconds, which increased the rate 
substantially above orbital rate and caused a slow drift in the opposite 
direction. Also in maneuver 3, an unexplained pitch-up, pitch-down 
perturbation occurred just before the maneuver was started. 
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Hi RAP Accelerometer Flight Data 

The HiRAP flight data were also recorded during the maneuvers and 
processed on the same postflight OEX-CCT as the rate gyro data. However, no 
postflight calibration of the HiRAP data was performed at the ODL, and the data 
received at LaRC were in the form of raw counts from the HiRAP data 
electronics. The counts from the X- and Z-axis accelerometers were converted 
to units of 10“ 6 g and plotted (relative to the initial value) in figures 
5-10. A steep negative slope caused by a large temperature-dependent bias is 
clearly evident in the figures. The actual accelerations, both X- and Z- 
channels, are known to be small, on the order of 10 -6 g, during the 
pre-maneuver hold period. During the maneuver, the theoretical induced 
acceleration in the X-channel based on the measured pitch rate is about 43x10 
g for maneuver 1 and about 21xl0‘ 6 g for maneuvers 2 and 3. The corresponding 
values for the Z-channel are -28xl0" 6 g and -13xl0 -6 g. These induced 
accelerations can be seen in the plots even in the presence of the temperature 
slope and large noise levels. 

Rotational! v Induced Acceleration 

A linear accelerometer mounted on a translating, rotating body at some 
distance from the body CG will have an output proportional to the algebraic 
sum of the projections onto the accelerometer axis of the linear acceleration 
of the body, the rotational ly Induced acceleration of the proof mass, and the 
gravity-gradient-induced acceleration of the proof mass. Previous analysis of 
the HiRAP data 5 Involved the extraction of the linear contributions which 
required an analytical correction for rotation based on the ACIP rate gyro 



measurements of rotation rates. The equations for rotationally induced 
acceleration 6 for the three HiRAP axes are 

a = -(q 2 + r2) X x + (pq - f) Y x + (pr + q) \ 
r x 

a = (pq + r) X + (p 2 + r2) Y y + (qr - p) Z y 

y 

a f = (pr - q) X z + (qr + p) Y z + (p 2 + q 2 ) T- z 

where a Pi , i = x, y, z, are the total rotationally induced accelerations; 
p, q, and r are angular velocities about X^, Y^> and Z^; and X-j , Y^ , 

Zj are offsets of the X, Y, and Z accelerometers from the CG along the ith 
body axis. The offsets for the X- and Z-channel accelerometers are shown two 

dimensionally (i = x, z) in figure 11. 

An in-depth analysis of the relative contribution of the various terms, 
based upon STS 61-C attitude maneuver data, was performed to aid simplification 
of these equations. The Y-axis was omitted from consideration because the 
orbiter CG and the Y-channel accelerometer are both very close to the lateral 
plane of symmetry so that the anticipated uncertainty in determination of the 
offset is of the same order of magnitude as the offset. The ACIP on STS 61-C 
had an inoperative p-channel rate gyro; however, average values of p were 
determined from the postflight attitude and trajectory history 7 for STS 61-C. 
These values are shown for comparison with the average ACIP q and r values in 
tables 1 and 2 for the pre-maneuver and maneuver segments, respectively. 
Retaining only the terms which produce accelerations greater than 0.1x10 g 

and assumming that p = q = r = 0 reduces the equations to 

a r - - (q 2 + r 2 ) \ 



8 


Gravity Gradient Effects 

The substantial distance (3.7 m) between the CG and the HiRAP mounting 
location contributes to non-negl igible gravity-gradient-induced acceleration in 
both X- and Z-channels. This acceleration is given for a circular orbit (valid 
assumption for STS 61-C) by 

1 ^ 1 

a = 2 w 2 (X' 2 + Z' 2 ) 1/2 | cos 0 cos (6 + tan — ) | 

q x z vi 

x *x 


-1 z ' 

a . (X' + Z; 2 ) 1/z | sin 6 cos (8 + tan t 2 ) I 

9 z 


where w is the orbital angular velocity (■ 0.0658 deg/sec for STS 61-C), 8 is 

I _ I 

the angle between the X& axis and local vertical, and 5( x and Z z refer 
to pre-flight values. The maximum value of ag x is 0.91xl0- 6 g, which occurs 
at 0 = -16.7° and 163.3°, and the maximum magnitude of a gz occurs at 

0 = -61.7° and 118.3° with a Qz = -.7xl0~ 6 g. Since each a Qi is a small 
contribution (5 percent or less) to the total measured acceleration, the use of 
x‘ and l z introduces negligible error. 

The gravity gradient torque exerted on the orbiter during the OEX 
maneuvers was primarily about the Yb axis and is given by 


V 


3 P e 

(I - I ) sin 0 cos 0 

r 3 xx zz 


The induced angular acceleration is then 


q q = I 


^ - 3 ■ 


yy 


R 3 i 


(I - I ) sin 0 cos 0 
v xx zz' 


yy 


9 


which reduces to approximately 

q = -3.6x10”® sin 0 cos 6 (sec -2 ) 

g 

for the STS 61-C orbiter mass properties and orbital altitude. This angular 
acceleration tends to stabilize the orbiter axis along the local vertical, 

o 

either nose up or nose down. The maximum absolute value of q g , 

1.8xl0” 6 sec” 2 , occurs at 9 = ±45° and ±135° and induces an X-channel linear 

acceleration of about ±0.36x10”® g and a Z-channel value of about 

±0.56x10”® g. The presence of this effect Is clearly evident in figure 12, 

which shows the maneuver 3 ACIP pitch rate history from 0 = 0° to 180° compared 
with the theoretical pitch rate obtained by integrating the equation for q g . 

Aerodynamic Forces and Moments 

The STS 61-C mission was flown during minimum solar activity when the 

atmospheric density was minimum at the 61-C orbital altitude. The expected 

aerodynamic acceleration for this condition, based on the latest orbiter 

flight-derived free-molecule-flow aerodynamic coefficients 8 , is less than 

0.1x10"® g and is considered to be negligible. The maximum aerodynamic 

pitching moment (from C. Cooke and others at Charles Stark Draper Laboratory, 

study of Shuttle experiment acquisition, tracking and pointing for Jet 
Propulsion Laboratory) produces a q of 2.3xl0" 7 sec -2 , which would induce 

negligible linear accelerations of 0.05x10"® g in the X-channel and 0.07x10 ® g 
in the Z-channel . 


HiRAP Zero Offset and Temperature Bias 
The outputs of the HiRAP accelerometer are temperature dependent and must 
be calibrated accordingly for each data set to provide the desired absolute 
measurement accuracy. This calibration is normally accomplished through 
analysis of flight data from an on-orbit calibration sequence.® The zero bias 
is established for a reference point, t re f, and the temperature bias slope in 
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integrated over the subsequent measurement interval based on measurements from 
accelerometer-mounted thermocouples. For a linear variation in temperature 
with time, this correction is equal to 

*T x -\ + B x - ‘ref' + C x<‘ - ‘ref’" 

a T z - A z + B z (t ' W * C z (t - l ref> 2 

which is the form used in this analysis. However, no separate calibration 
analysis was performed to solve for the coefficients A-j , B-j , and C-j ; 
rather, their solution was obtained simultaneously with the solution for the CG 
offsets. 


Solution for 


From the foregoing considerations, 
the significant sources are 

a - - (q 2 ♦ r*) X x ♦ a ♦ 5 g 2, 

X A 


CG Offsets 

the total induced accelerations from 


The ith accelerometer flight measurement is then equal to these induced 
accelerations plus the reference zero offset and temperature bias at time t. 
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These equations were used with a least squares routine with matrix Inversion to 
solve for the four unknowns in each equation; X x , A x , B x , and C x for 
the the X-channel and Z z , A z , B z , and C z for the Z-channel; while 
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minimizing the difference between the measured acceleration on the left and the 
calculated acceleration on the right. The data sample used in this solution 
includes a period during the pre-maneuver attitude hold and an interval during 
the maneuver. The solution is equivalent to solving for the CG offsets based 
on the change in induced acceleration for a given change in the pitch angular 
velocity between the two periods, an approach which is forced by the lack of 
absolute calibration of the measured accelerations. 

Selection of Data Segments for Analysis 
Figure 13 shows graphically the definition of the two data segments used 
for each maneuver in the CG offset calculations. The pre-maneuver hold segment 
from tj to t 2 was selected to be free from major unmodeled perturbations 
and to end prior to the thruster firing which initiated the maneuvers. The 
maneuver segment from t 3 to was started after the thruster firing 
accelerations had damped to a negligible level and contained no major unmodeled 
perturbations. The selected segments for the three maneuvers are listed in 
table 3 with the values for t re f which were taken to be the approximate time 
of maneuver initiation. Each of these segments is approximately 100 seconds 
long and contains about 11,200 measurements for a total of about 22,400 for 
each maneuver. 


CG Offset Solutions 

Solutions were obtained for X x and Z z for the three maneuvers using 
the data segments defined in table 3 and also for a group of restricted data 
sets. The restricted data sets were defined by applying the criterion that any 
measurement whose residual (measured acceleration minus calculated 
acceleration) in the unrestricted solution exceeds no is deleted from the data 
set. A value of n greater than 7 produces essentially an unrestricted data 
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set, while decreasing values of n produce data sets which are increasingly more 
restricted to the central measurements. The CG offsets obtained for the 
unrestricted data and for the sets with n = 6, 5, 4, 3, and 2 are listed in 
table 4 and also shown in figure 14 compared with the pre-flight values. 

Significant variations in X x and Z z are present in the solutions from 
maneuver to maneuver and also across the restricted solutions for a given 
maneuver. These results show a strong sensitivity to the distribution of 
"noise" in the data, which indicates violation of the assumption of random 
normal distribution of noise. The distributions of residuals over the 
pre-maneuver and maneuver segments are shown for the unrestricted solutions in 
figures 15-20 with the theoretical normal distribution curve shown for 
comparison. Some significant departures from random normal are seen in these 
distributions, both in the central values which would account for the variation 
in the CG offset solution from maneuver to maneuver and in the outlying values 
which would account for the variations in the restricted solutions. Without 
exception, the quality of fit of the flight data obtained in each solution 
varies between the two segments with the most pronounced difference occurring 
in the maneuver 3 X-axis solution, figure 17, for which the far better fit of 
the maneuver segment is evident in the high concentration of central values. 

In general, these distributions show much larger than expected (based on the 
random normal assumption) concentrations of outlying values above 3o which are 
typically non-symmetrical about the mean. Further, the outlying values are not 
randomly distributed over the data segment but occur as one or more large 
peaks. For example, the maneuver 3 Z-axis pre-maneuver segment contains 43 
measurements whose residuals exceed 4o, only 3 of which are negative, with 9 
positive values in one single peak and 30 positive values In a second peak. A 
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similar situation exists with the maneuver segment for this case which has 55 
measurements above 4o, of which 46 are positive and contained within four 
peaks. The probability of occurrence of values above 4o for a random normal 
distribution is only 6xl0 - ^, or one value in the 22,400 measurements compared 
with a total of 98 for this case. 

Analysis of Noise/Vibration Effects 
The noise level on the rate gyro data is about 4 to 6 percent of the 
measured signal and appears to be random. However, the noise level in the 
accelerometer data is of the same order as the rotational ly induced 
acceleration and frequently exceeds ±50xl0 -6 g. Figure 21 shows the X-channel 
accelerometer output for an expanded time scale which depicts the variations 
more clearly than does figure 5. The signal does not exhibit the 
characteristics of random noise but rather reflects the expected orbiter 
structural vibration characteristics; e.g., a composite waveform produced by 
several vibration sources such as motors, pumps, and directional antenna 
systems. 

A fast Fourier transform routine was used to determine the power spectral 
density in the accelerometer output for each of the data segments as shown in 
figures 22-27. The sample size for this routine was 8192 (2^), reduced from 
the segment size of 11,200, and therefore covers only 73 seconds of the 100- 
second segments. Noticeable differences between the spectra are seen when 
comparing the pre-maneuver and maneuver segments for maneuvers 1 and 3 in both 
X- and Z-channels. For example, the pre-maneuver sequence during maneuver 3 
for the X-channel (figure 24(a)) shows much lower amplitudes below 2 Hz than 
the maneuver segment in figure 24(b), while the opposite is true of maneuver 1 
in figures 22(a) and (b). More dramatic differences are seen in the Z-channel 

1 and 3 where large amplitude spikes are seen in one 


spectra for maneuvers 
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segment (figures 25(a) and 27(b)) but not in the other segment. Also, the 
orbiter structural resonance near 5 Hz (first longitudinal bending mode) is 
stronger in the Z-channel. Given the changes that frequently occur between the 
pre-maneuver and maneuver segments, significant changes could be expected even 
during the span of a single segment. Such changes could certainly alter the 
degree to which the vibration signal is randomized and thereby affect the CG 
offset solution. A model of the vibration environment for the maneuver 3, 
Z-axis case was developed from its component spectrum and was used to generate 
a simulated vibration signal which was superimposed on the theoretical output 
of HiRAP for the target orbiter rotation rates. Sine waves spaced at 
0.0138 Hz from 0.1 to 3.0 Hz with amplitudes consistent with the component 
spectra were started at zero phase angle 50,000 sec prior to t re f in order to 
randomize the composite signal. The simulated HiRAP data were input to the CG 
offset solution routine which produced the residual distribution shown in 
figure 28. This distribution shows more symmetry and a higher concentration 
near the mean when compared with the distributions from the flight data 
solutions. 


Unfiltered Orbiter Vibration 

The power spectral density profiles obtained from the HiRAP data contain 
few clues, if any, to the nature and source of the vibrations which would lead 
to improvement of the CG solutions. The possibility that the low-frequency 
0 to 3 Hz spectrum could be in part due to beating of closely spaced 
frequencies above 5 Hz led to the generation of power spectral density profiles 
which were corrected for the 43 db per decade rolloff of HiRAP between 2 and 20 
Hz. These spectra are shown in figures 29-34. The first longitudinal bending 
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mode resonance near 5 Hz is more pronounced and the large directional antenna 
"dither" motion is seen at 17 Hz. Higher frequency components in the range of 
5 to 15 Hz are much stronger in the Z-channel as compared with the X-channel. 

No obvious correlation between the higher frequencies and the lower frequencies 
is present. These data indicate, within the limits of the fast Fourier 
transform technique^, the expected level of vibrational ly induced acceleration 
at the HiRAP mounting location, even during a so-called "quiet" period, which 
would bear on any "zero-g" investigation or process. 

Conclusions 

The technique for determination of the orbiter center of gravity from 
flight measurements as described herein produces results which are reasonably 
accurate considering the complex acceleration environment in which the 
measurements were obtained. The orbiter structural vibration produces a 
complex acceleration waveform which is composed of a number of different 
frequencies generated by various mechanical subsystems. The operating schedule 
for the subsystems varies as does the degree to which the resultant waveform is 
randomized. The occasional in-phase amplitude peaks are not random over the 
relatively short time segments used in this analysis. Consequently, the 
center-of-gravity solution is sensitive to the choice of data segment length 
and position and also to any statistical culling procedure. The observed 
variation in the solutions is due primarily to variations in the sample 
distribution and not to measurement errors or first-order errors in the 
temperature bias model; therefore, any improvement in measurement accuracy or 
temperature model would not improve the solution. Such improvement would 
require more sophisticated analytical and/or statistical methods beyond the 
scope of this study. The improved low-pass filter on the OARE will further 
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attenuate the vibration signals and the digital filtering techniques to be 
applied during the processing may further reduce the vibration-induced errors 
in the center-of-gravity solution. 
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Table 1 Average rotational rates during pre-maneuver hold (deg/sec) 



Maneuver 1 

Maneuver 2 

Maneuver 3 

p 

.001 

-.0013 

.002 

q 

.06 

.06 

.1 

r 

-.014 

-.013 

-.011 



Table 2 Average rotational 

rates during maneuvers (deg/sec) 


Maneuver 1 

Maneuver 2 

Maneuver 3 

P 

-.01 

.018 

.008 

q 

-.67 

-.47 

-.47 

r 

-.029 

.0036 

-.0014 


Table 3 Data Segments (GMT seconds) 



Maneuver 1 

Maneuver 2 

Maneuver 3 

l l 

31560 

60068 

55675 

l 2 

31660 

60168 

55775 


31700 

60177 

55805 

t REF 

31730 

60210 

55840 


31830 


60310 


55940 
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Table 4 CG offsets obtained for various data restrictions. 


X x (pre-flight value = -3.07 m) 

Maneuver 1, Maneuver 2, Maneuver 3, Mean 

Restrictions m m m ± deviation, m 



none 

-2.932 

-3.112 

-3.077 

-3.040 ± .078 


6o 

-2.956 

-3.132 

-3.205 

-3.098 ± .105 


5o 

-2.968 

-3.132 

-3.285 

-3.128 ± .129 


4o 

-2.982 

-3.123 

-3.277 

-3.127 ± .120 


3o 

-3.019 

-3.118 

-3.103 

-3.080 ± .044 


2o 

-3.142 

-2.987 

-3.063 

-3.064 ± .063 




(pre-flight value 

= 2.021 m) 



Restrictions 

Maneuver 

m 

1, Maneuver 2, 

m 

Maneuver 3, 
m 

Mean 

± deviation, m 


none 

1.986 

1.894 

2.145 

2.008 ± .104 


6o 

2.058 

1.971 

2.145 

2.058 ± .071 


5o 

2.030 

1.956 

2.107 

2.031 ± .062 


4o 

2.004 

1.942 

2.070 

2.005 ± .052 


3o 

1.953 

1.933 

2.125 

. 2.004 ± .086 


2o 

1.974 

1.931 

2.118 

2.008 ± .080 

- 



- - 









1 a i ned 





Figure 5. H i RHP X-channel relative output, maneuver 




Figure G. HiRflP X-channel relative output, maneuver 




Figure 7. HiRHP X-channel relative output, maneuver 




Figure 8. HiRRP Z-channe 1 relative output, maneuver 





Figure 9. HiRHP Z-channel relative output, maneuver 




HlRRP Z-channel relative output, maneuver 



orb i ter 



Figure 11. Convention -for measurement of CG offsets. 
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Figure 12. Effect of gravity-gradient-induced acceleration on 
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Figure 13. Description of data Intervals for CG solution. 
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H1RRP location with respect to CG, 
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(a) Pre-maneuver segment 

Figure 15. Distribution of residuals, maneuver 1, 
X-axis unrestricted solution. 
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(b) Maneuver segment 
Figure 15. Concluded. 
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(a) Pre-maneuver segment 

Figure 16. Distribution of residuals, maneuver 2, 
X-axls unrestricted solution. 
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(b) Maneuver segment 
Figure 16. Concluded. 



F r act i on 
of total 
measurements 



(a) Pre-maneuver segment 

Figure 17. Distribution of residuals, maneuver 3, 
X-axis unrestricted solution. 
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Figure 17. Concluded. 
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(a) Pre-maneuver segment 

Figure 18. Distribution of residuals, maneuver 1, 
Z— axis unrestricted solution. 
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Figure 18. Concluded. 
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(a) Pre-maneuver segment 

Figure 19. Distribution of residuals, maneuver 2, 
Z-axts unrestricted solution. 
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(b) Maneuver segment 
Figure 19. Concluded. 
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Ca) Pre-maneuver segment 

Figure 20. Distribution of residuals, maneuver 3, 
Z— axis unrestricted solution. 


Fr act 1 on 
of total 
me asu rements 



Figure 20. Concluded. 




H1RRP X-channel relative output, maneuver 1, expanded scale. 
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(a) Pre-maneuver segment 
Figure 22. Power spectral density, maneuver 1, X-axis 
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(b) Maneuver segment 
Figure 22. Concluded. 
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Figure 23. 
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Cb) Maneuver segment 
Figure 23. Concluded. 
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(a) Pre-maneuver segment 
24. Power spectral density, maneuver 3, X-axis 



(b) Maneuver segment 
Figure 24, Concluded. 
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(a) Pre-maneuver segment 

Figure 25. Power spectral density, maneuver 1, Z— axis. 
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(b) Maneuver segment 
Figure 25. Concluded. 
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Figure 26. 
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Power spectral 


Pre-maneuver segment 
density, maneuver 2, Z— axis. 
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Cb) Maneuver segment 
Figure 26. Concluded. 
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(a) Pre-maneuver segment 

Figure 27. Power spectral density, maneuver 3, Z-axfs 
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(b) Maneuver segment 
Figure 27. Concluded. 
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(a) Pre-maneuver segment 

Figure 28. Distribution of residuals for simulated vibration 
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(b) Maneuver segment 


Figure 28. Concluded. 
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Figure 29. 
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(a) Pro-maneuver segment 

Figure 30. Un-f titered power spectral density, maneuver 2, X-axis. 
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(b) Maneuver segment 
Figure 30. Concluded. 
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(a) Pre-maneuver segment 

Figure 32. Unflltered power spectral density, maneuver 1, Z-axIs 
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(b) Maneuver segment 
Figure 32. Concluded. 
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(a) Pre-maneuver segment 

Figure 34. Unfiltered power spectral density, maneuver 3, Z-axis 


-7 



Frequency, Hz 


(b) Maneuver segment 
Figure 34. Concluded. 



2. Government Accession No 


NASA Report Documentation Page 

•V'OTkYjV •• HfXl ' 

st ,Kv Vr »' •$*' tfmf 

1 Report No^ " 2. Government Accession No 

NASA TM-102627 

4. Title and Subtitle 

Determination of Shuttle Orbiter Center of Gravity 
from Flight Measurements 


3. Recipient's Catalog No. 

5 Report Date 

January 1991 

6. Performing Organization Code 


7. Author(s) 

E. W. Hinson, J, Y. Nicholson, and R. C. Blanchard 


9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665-5225 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


8. Performing Organization Report No. 

10. Work Unit No. 

506-40-91-01 

11. Contract or Grant No. 

13 Type of Report and Period Covered 

Techni cal Memorandum 

14 Sponsoring Agency Code 


15. Supptementary Notes 

E. W. Hinson: ST Systems Corporation, Hampton, Virginia. 

J. Y. Nicholson: Vigyan Research Associates, Inc., Hampton, Virginia. 

R. C. Blanchard: Langley Research Center, Hampton, Virginia. 


16. Abstract 

Flight measurements of pitch, yaw, and roll rates and the resultant rotationally 
induced linear accelerations during three orbital maneuvers on Shuttle mission STS 
61 -C have been used to calculate the actual orbiter center-of-gravity location. The 
calculation technique reduces error due to lack of absolute calibration of the 
accelerometer measurements and compensates for accelerometer temperature bias 
and for the effects of gravity gradient. Accuracy of the technique was found to be 
limited by the nonrandom and asymmetrical distribution of orbiter structural vibration 
at the accelerometer mounting location. Fourier analysis of the vibration was 
performed to obtain the power spectral density profiles which show magnitudes in 
excess of 10^ ug^/Hz for the actual vibration and over 500 ug^/Ftz Tor the filtered 
accelerometer measurements. The data from this analysis provide a 
characterization of the Shuttle acceleration environment which may be useful in 
future studies related to accelerometer system application and "zero-g” 
investigations or processes. ^ 


17 Key Words {Suggested by Author(s)) 

^ — — 

18. Distribution Statement 


Accelerometer 


Unclassified-Unlimited 


Center-of-gravity 



- 


Rotational dynamics 


Subject 

Category 34 


19 Security Classif. (of this report) 
Unclassified 

20 Security Classif. (of this page! 

Unclassif ied 

21 No. of pages 

54 

- 

22. Price 

A0 4 


NASA FORM 1826 OCT 86 



